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Abstract 

We calculate the I-V characteristics of a tunnel junction containing impurities 

in the barrier. We consider the indirect resonant tunneling involving the 

impurities. The Coulomb repulsion energy E c between two electrons with 

opposite spins simultaneously residing on the impurity is introduced by an 

Anderson Hamiltonian. At low temperatures T <C E c the I-V characteristic 

is linear in V both for V < E c and for V > E c and changes slope at V = E c . 

This behavior reflects the energy spectrum of the impurity electrons - the 

finite value of the charging energy E c . At T ~ E c the junction reveals an 

ohmic-like behavior as a result of the smearing out of the charging effects by 

the thermal fluctuations. 

Ms. No PACS numbers: 73.20.Dx, 73.40. Gk. 
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In this paper, we study the transport through a tunnel junction with a thin, amorphous 
insulating layer separating the metal leads. At low temperatures two tunneling channels 
contribute to the conductance. The first one is the direct tunneling between the leads. 
The second one is a resonant tunneling through impurity states present in the barrier. This 
mechanism is dominant for not too small thickness of the layer For higher temperatures 
and larger applied voltages the inelastic tunneling involving pairs of localized states with 
phonon emission or absorption becomes important [|3|-||. We restrain our study to the case 
of resonant tunneling and calculate the current-voltage characteristics of such junctions. 
The current flows through the junction by tunneling of electrons from the emitter electrode 
to the impurity and then to the collector. We calculate the tunneling current through the 
impurity by taking into account the Coulomb repulsion energy between two electrons with 
opposite spins simultaneously residing on the localized state. In systems such as a-Si tunnel 
junctions this energy is rather large (0.1 — 0.2 eV) [||||. Experimentally measured current 
for a given dc voltage bias is obtained by averaging the currents through the individual 
impurities over their positions and energy distribution. 

The Hamiltonian of the system under consideration consists of three terms: 

Htot — Pleads + H imp + H tun . (1) 

The Hamiltonian describing the leads is written as: 

Hl eads = 6 k a L a ka + e p h l° h P°- ( 2 ) 
ka po- 

The impurity Hamiltonian is: 

Himp = e c c i c ° + E c n ^i, (3) 

and the tunneling term is: 

Htun = ^( T ^4 a fccr + H.c.) + Y( T Rp h la c <r + H.c). (4) 

ka Pcr 

The Hamiltonian is expressed in terms of annihilation(creation) operators for the emitter, 
b prT for the collector and c a for the impurity with k(p) the corresponding quasimomenta and 
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o =|, I is the spin index. e£ and are the single-particle energies in the emitter and the 



collector leads, respectively. E c is the Coulomb repulsion between electrons with opposite 
spins. The single-electron energies are measured from the corresponding Fermi levels \xl 
and fiR in the emitter and collector, and the dc bias is fiL — = eV. The resonant 
level energy is e c = e° c + aeV, where e° c is the bare resonant level energy and a measures 
the portion of the voltage drop on the localized center. The particle-number operator is 
ricr = c\c a . In Eqn. (gj) for the tunneling Hamiltonian, Tlu and Tr p are the tunneling matrix 
elements. Note that they depend on the impurity position in the junction. We assume that 
the impurities are randomly distributed in the junction barrier and their resonant levels are 
uniformly distributed in energy. 

First we calculate the impurity's Green's function. To correctly account for the non- 
equilibrium nature of the system we apply the Keldysh technique |7j . In this technique one 
introduces the retarded (advanced) and distribution Green's functions. We further assume 
that the relaxation processes in the leads are much faster than in the impurity. Thus, we 
can consider the leads as equilibrium systems and the corresponding Green's functions are 
given by the usual expressions for a non-interacting equilibrium electron systems 0. Here 
we give only the final expressions for the impurity's Green's functions (for details on the 
derivation see f|). The retarded Green's function is obtained in the form 



to — e ca — £ — £1 



(5) 



(oj - e c - S ) (a> - e c — £ - £i) - E c {u - e c - E - £2) ' 



where the analytical continuation ui 



oj + i0 + is used. The self-energy parts in Eqn. ([|) 



are given by 



k,i=L,R 




(6) 
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where f{e\) is the Fermi-Dirac distribution function for the electrons in the leads. We have 
used the following notations: e ccr = e c + E c < n_ a > and e ca = e c + E C {1— < >). This 
Green's functions describes two levels for the impurity electrons - a lower, resonant level 
with energy e c and the upper level with energy e c + E c . 

The distribution Green's function is calculated under the assumption that the transient 
processes after the switching-on the dc bias have decayed and the result is 



E|T Lfc | 2 A < (A;,^) + E|^ p | 2 J B < (p, 



<[UJ> = E \T Lk \\A r {k,uu) - A a (k,u)) + E \T Rp \ 2 (B r (p,u) - B a (p,u)) 

k P 

(G r (u) - G a {uj)) , (9) 

where A r ^ a>< (B r ^ a>< ) are the retarded, advanced, and distribution Green's functions in the 
left (right) lead. 

The average number of impurity electrons is calculated by solving the following equation 

n = -J^-ImG < (u) (10) 

(we assume < n a >=< n_ a >= n). 

It should be mentioned ( and references therein) that the Green's function formalism 
presented here is valid only for temperatures higher than the characteristic temperature for 
this problem - the Kondo temperature Tk- 

Next we calculate the current through the impurity. In a steady-state it is given by 

J 2tc >y L (uj) + >y R {uj) 
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where 



lL (co) = 7rJ2\T L k\ 2 6(u J -e L k ), lR {u) = n £ \T Rp \H{u - ef) (12) 
fc p 

are the elastic widths of the impurity level due to tunneling through the left (right) barrier, 

respectively. We take broad, flat density of states for the leads' electrons and in the following 

7l(c<j) and 7r(co) are assumed to be constants, 7^ and 7^, independent of energy. 

Now we compare our treatment with the work of Glazman and Matveev || . Their results 
can be obtained from ours by approximating the well electrons Green's function (Eqn. (gD) 
with the corresponding expression for an isolated one-center Hubbard model |T(| i.e. one 
has to set all the self-energies £i(co>), i = 0,2 (Eqns. (6-8)) to zero. This approach neglects 
the virtual tunnel processes and it is justifiable when the temperature T is much larger 
than the elastic level width 7 = 7l + 7r- With this approximation the integral equation for 
the average number of well electrons (Eqn. 10) reduces to an algebraic one whose solution 
coincides with Eqn. 6 in Ref. 6. Similarly, formula (11) for the current gives the simple 
analytical expression Eqn. 8 in Ref. 6. 

To obtain the current through the junction we must integrate over the position of the 
impurity and to average over the energies of the impurity levels. If we denote the position 
of the impurity, relative to the left lead, by x then 7^ and 7^ depend exponentially on x 

( 2(d-x)\ 
7i? = 7oexpl I, (13) 

where d is the width of the tunnel junction and £ is the impurity localization length. In a 
first approximation the coefficient 70 can be considered independent on x. One can show 
that most significant contribution to the current through the junction give the impurities 
which are close to its middle - x = d/2. Therefore, we simplify the calculation by taking 
x = d/2 in Eqn. (|H|) and discarding the impurities which are far from this optimal position. 
With this approximation the current through the junction is given by 

I = SdJ de c g(e c )I(e c ) (14) 



7 L = 70 exp 
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where S is the area of the junction and g(e c ) is the density of impurity states. In further 
calculations we take g(e c ) to be a constant. In Eqn. ( |14D the integration is over all impurities 
which give non-zero contribution to the current. 

In Fig. 1 we show the current through the impurities as a function of the energy of the 
impurity level for: a) V = 2.2E C ; b) V = 0.6E C and for T = 0.08£ c (70 = 0.02£ c ). We have 
set fin = since the results are independent on the value of the right chemical potential. 
The current is almost zero when both impurity levels - the resonant and the upper level, 
are well below the right chemical potential or well above the left chemical potential. In the 
former case n — 1 - there are 2 electrons on the impurity and in the latter n ~ 0. For V > E c 
the current is maximal when both impurity levels are between the left and right chemical 
potentials - there is approximately 1 electron on the impurity. The left "shoulder" in I(e c ) 
dependence for V = 2.2E C (Fig. 1 a) corresponds to impurities for which n « 1/3 and the 
right "shoulder" is due to impurities with n « 2/3. For V < E c there are two peaks in the 
J(e c ) dependence (Fig. 1 b). The lower comes from tunneling through the upper level for 
impurities deep below jj, R (for them only the upper level is above /ir). The upper peak is 
due to impurities whose resonant level is above /ir and the upper level is above //£,. In this 
case the tunneling current flows predominantly through the resonant level. In Fig. 1 c we 
show I(e c ) dependence for T = 0.9E C and V = 0.6E C . It is evident that for temperatures 
comparable with E c all charging effects are smeared out by the thermal fluctuations. 

In Fig. 2 we present the I-V characteristic of the tunnel junction, calculated for: a) 
T = 0.002£' c and b) T = 0.9E C . In the former case the current is linear in V both for low 
voltage V < E c and for V > E c and changes slope at V = E c - the slope diminishes for 
V > E c . This behaviour can qualitatively be understood by the following considerations. 
With increasing V more impurities are involved in the transport through the junction. As V 
increases to E c the current through the junction increases significantly since it flows through 
impurities which give maximal current (Fig. 1). Increasing V above E c includes impurities 
with resonant levels close to - for them the current is rapidly decreasing function of 
e c (Fig. 1). We must point out that for temperatures T < 0.5E C the I-V characteristic 
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is practically independent on the temperature which again is a signature of the charging 
effects. 

The I-V characteristic at T = 0.9E C (Fig. 2b) shows ohmic behavior of the tunnel 
junction. The current is linear in V for every voltage - the charging effects of the resonant 
tunneling are fully masked by the thermal fluctuations. 

In this work we have approximated the tunneling matrix elements to be independent on 
energy (equivalently, on the applied dc voltage). One may speculate that taking T L ^ to 
be increasing functions of V and e c (softening of the barrier with increasing energy) will 
give the experimentally observed I-V characteristics H - for V exceeding E c the slope of the 
/(V) dependence increases. 

In summary, in the present paper the phenomenon of resonant tunneling through a tunnel 
junction containing impurity states in the barrier region has been studied with the Coulomb 
repulsion between electrons on the impurity taken into account. The I-V characteristic at 
temperature T <C E c is linear in V both for V < E c and for V > E c and its slope diminishes 
at V = E c . The tunnel junction has ohmic-like I-V characteristic for temperatures of 
the order of the charging energy when the charging effect is smeared out by the thermal 
fluctuations. 
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FIGURES 

FIG. 1. Current (in arb. units) through the impurity centers as a function of their energy e c 
for T = 0.08£ c : a) V = 2.2E C (the dotted line), b) V = 0.6E C (the solid line), and c) for T = 0.9£ c 
and V = 0.6E C (the dashed line). 

FIG. 2. I-V characteristics of the tunnel junction for: a) T = 0.002£' c (the solid line), and b) 
T = 0.9E C (the dashed line). The current is normalized to I(E C ). 
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